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We demonstrate a versatile method to create state-dependent optical lattices by applying a mag¬ 
netic field gradient modulated in time. This allows for tuning the relative amplitude and sign of 
the tunnelling for different internal states. We observe substantially different momentum distribu¬ 
tions depending on the spin-state of fermionic 40 K atoms. Using dipole-oscillations we probe the 
spin-dependent band structure and find good agreement with theory. In-situ expansion-dynamics 
demonstrate that one state can be completely localized whilst others remain itinerant. A systematic 
study shows negligible heating and lifetimes of several seconds in the Hubbard regime. 
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Ultracold atoms in optical lattices provide a highly 
tunable platform to simulate the behaviour of electrons 
in solids. When the tunnelling in the lattice depends 
on the internal spin state, SU(2) symmetry is explicitly 
broken and novel quantum phases emerge, such as un¬ 
conventional superconductivity owing to a Fermi-surface 
mismatch, or exotic forms of magnetism arising from 
anisotropic spin exchange [1, 2]. Both the static proper¬ 
ties and the dynamics of particles in the lattice will then 
depend on their internal state. Realizing such a state- 
dependent tunnelling in an optical lattice requires a cou¬ 
pling between internal and external degrees of freedom 
and has so far been demonstrated with bosonic atoms, 
allowing for quantum computation and simulation [3-13]. 
The technique used there mainly relies on the differen¬ 
tial coupling of the lattice laser field to different atomic 
transition lines [14]. Its range of applicability is however 
limited by the intrinsic problem of heating by sponta¬ 
neous emission. This problem is particularly severe for 
fermionic potassium and lithium atoms because of their 
small fine-structure splittings, which has hindered the re¬ 
alisation of state-dependent optical lattices for fermions. 
Other proposed methods, relying on earth-alkaline and 
similar atoms [15, 16], or on atom chips [17], seem to in¬ 
volve difficult experimental challenges. 


Here we present the implementation of a spin- 
dependent lattice for ultracold fermions using a different 
method. Following the proposal in Ref. [18], our method 
relies on the application of an oscillating force to the 
particles in the lattice, with an amplitude which depends 
on their internal state [19]. The resulting system can 
then be well described by an effective time-independent 
Hamiltonian [20-23] , with differently renormalized tun¬ 
nelling terms for each internal state. The general idea 
can readily be extended to mixtures of different atomic 
species or other artificial or conventional lattices. In our 
case the spin-dependent force in a particular direction x 
on state <r is provided by a magnetic field gradient, and 



FIG. 1. (a) Magnetic moments of the F = 9/2 hyperfine 

manifold of 40 K as a function of the external field. Blue, red 
and green denote the \rriF = —9/2), |—5/2) and | — 1/2) sub- 
levels throughout, (b) The effective energy bands in quasimo¬ 
mentum space are given by the time-average of the “shaken” 
bands. As the force depends on the internal state, so does 
the effective band, (c) A cloud of 40 K (blue) is trapped in 
a retro-reflected laser beam. An oscillating current I(r) in a 
single coil creates the oscillating gradient cU|B|. The uniform 
external field B cx t is provided by additional coils. 


is given by 

where E a gives the energy of a state as a function of the 
magnitude of the external magnetic field |B| [24]. For 
the F = 9/2 hyperfine-manifold of 40 K used in our ex¬ 
periment, the force resulting from a given gradient can 
take on various positive or negative values or vanish, de¬ 
pending on the Zeeman-sublevel, see Fig. la. 

An oscillating force renormalizes the amplitude and 
phase of tunnel-couplings in a lattice, which has been ob¬ 
served in optical lattices for ultracold bosons [20, 25-28]. 
In a single-band tight-binding model, the time-dependent 
Hamiltonian of a one-dimensional non-interacting system 
is given by: 

H(t) = -t^c] ia c j+lta +H..c.-a'^2FZ(T)Y^3c],aC j ,v 
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where Cj a and Cj a are the creation and annihilation oper¬ 
ators of one fermion with spin a on site j. t is time, a the 
lattice constant and t the tunnelling energy. Using Flo- 
quet theory, an effective time-independent Hamiltonian 
can be derived, which describes the system on time-scales 
longer than the oscillation period [21, 22, 29-31]. For a si¬ 
nusoidally oscillating force F x [r) = Kg Hcos/a • sin (wsr), 
with angular frequency cos and a dimensionless mod¬ 
ulation amplitude Kg , the effective Hamiltonian reads 
(see [20-23] for a derivation): 

H eS = -tJ2 Jo(Kg) b- + Hc - (2) 
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The tunnelling has therefore been renormalized to a spin- 
dependent value 

Ks = tJo(Kg) (3) 

given by a 0 th -order ordinary Bessel function J 0 . An 
intuitive picture of this effect can be gained by consid¬ 
ering the time-dependent band energy as a function of 
quasimomentum q x in a co-moving frame. As illustated 
in Fig. lb, the average of the “shaken” band energy over 
one period then gives the effective band: 

4f (Qx) = 2f cos (aq x - a Fg(r)dr^ (4) 

Going beyond the non-interacting single-band regime, 
these effective Hamiltonians can contain additional terms 
such as longer-range tunnelling, as recently observed [18, 
32, 33], which could then also be made state-dependent. 

In a first experiment we study how the changing band 
structure of the effective Hamiltonian affects the quasi- 
momentum distribution of two different spin states. For 
this measurement, we prepare a degenerate cloud of 
1.9(4) x 10 4 non-interacting spin-polarized fermionic 40 K 
atoms in an optical dipole trap operating at 826 nm and 
load them into the lowest band of a one-dimensional op¬ 
tical lattice with a lattice constant of a = 532 nm and a 
tunnelling rate of t/h = 174(9) Hz. For details of the 
preparation procedure and trapping parameters, see [34], 
We then ramp up an oscillating current with a frequency 
of 750 Hz within 100ms, which runs through a single coil 
mounted about 1 cm away from the atoms, see Fig. lc. 
Reaching the first zero of t^ 2 at K 0 9//2 ss 2.4 requires 
a gradient amplitude of about 24G/cm corresponding to 
a current amplitude of 6.4 A. 

Fig. 2a shows the resulting quasimomentum distribu¬ 
tions in the lattice, measured using a band-mapping tech¬ 
nique [34] . The distribution broadens when the modula¬ 
tion is increased, because the width of the lowest band 
decreases. Above the critical value of K 0 9 ^ 2 ss 2.4, a 
double-peak feature appears for atoms in the mp = —9/2 
Zeeman sub-level (henceforth denoted |—9/2), and sim¬ 
ilarly for other mp values). This occurs because t f 9 ^ 2 



FIG. 2. Fermions in spin-dependent bands, (a) Quasimo- 
mentum distribution in the lattice, summed over q y and q z , 
as a function of gradient modulation amplitude for the |—9/2) 
(left) and | —5/2) (right) state. The effective band structures 
for each spin state are shown as diagrams for = 0, 1.2, 

2.4 and 3.0. (b) Second moment of the distributions. Blue 
circles (red diamonds) denote |—9/2) (|—5/2) atoms). Data 
show mean ± s.d. of 5 measurements, (c) Dipole oscillations 
in the spin-dependent bands (illustrated in the inset) show 
how the tunnelling is renormalized for each spin, relative to 
the oscillations without modulation (which have frequency 
vo)- The oscillation frequency squared is proportional to the 
tunnelling and is determined from fits to at least 60 measure¬ 
ments of the time-dependent quasimomentum peak position. 
Error bars show the fit uncertainty and solid lines the Bessel- 
functions for each spin, calculated without free parameters. 
The star indicates |—9/2) atoms without a lattice, (d) Tun¬ 
nelling ratio of the two spins as a function of modulation 
amplitude. 


becomes negative and therefore the band has minima 
at quasimomenta q x = ±7r /a rather than at q x = 0. 
The situation is very different when using the |—5/2) 
state. At the offset-field of 57.53(1) G used here, Kg°^ 2 = 
0.636(2) xKg 9 ^ 2 for the same magnetic field gradient am¬ 
plitudes. Therefore, when the |—9/2) atoms experience a 
completely flat band, those in the |—5/2) state still tun¬ 
nel at a rate of 86(4) Hz. As can be seen in Fig. 2b, their 
spread increases only slightly. The two states experience 
very different effective band structures, which allows for 
creating a tunable Fermi-surface mismatch. 
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In order to directly measure the effective tunnelling, 
we perform dipole oscillations in the modulated optical 
lattice. The tunnelling is proportional to the oscillation 
frequency squared, which, in the limit of small oscilla¬ 
tions, can be understood as arising from a change in ef¬ 
fective mass [37]. The atoms are displaced by 0.31(4) 7r/a 
in ^-direction and allowed to evolve for up to 350 ms. 
The oscillation frequency is then extracted from the 
time-dependent peak position of the quasimomentum- 
distribution. The maximum oscillation frequency of 
8.4(3) Hz is much smaller than the modulation frequency 
of 750 Hz, meaning that the dynamics should be well de¬ 
scribed by H c g. As shown in Fig. 2c, the oscillations for 
|—9/2) atoms become slower when the modulation ampli¬ 
tude is increased, as expected from the Bessel functions in 
Eq. 3. The |—5/2) atoms, on the other hand, experience 
a weaker modulation force and their oscillation frequency 
therefore changes much less. The spin-dependent oscil¬ 
lation frequency shows that the atoms behave as though 
they had different masses in the lattice. The effect is ex¬ 
pected to vanish in the absence of a lattice, as a quadratic 
dispersion is not changed by an oscillating force. Indeed, 
Fig. 2c shows that we observe no reduction of the oscilla¬ 
tion frequency when applying an oscillating gradient in 
a harmonic trap. 

By tuning the modulation amplitude, the tunnelling 
ratio of the two spins can be set to any positive or 
negative value, as shown in Fig. 2d. We now focus on 
the case where K 0 9,/2 ss 2.4 and the |—9/2) atoms 
experience a completely flat band with zero tunneling. 
They are therefore pinned to the lattice, whilst atoms 
in other states remain itinerant, see Fig. 3a. This sit¬ 
uation can be observed in the in-situ expansion of the 
atomic cloud. For these measurements we work with sin¬ 
gle spin states at a uniform offset-field of 208.15(1) G, 
as typically used for experiments with interacting 40 K. 
Starting from a trap configuration where w x Vt z = 27r x 
(67.8(3), 60.4(4), 233.5(3)) Hz we suddenly switch off the 
confinement along the lattice direction and measure 
the width of the cloud as a function of time [38, 39]. 
Fig. 3b shows that no expansion of the cloud can be ob¬ 
served for the |—9/2) state, where K^ 9 ^ 2 = 2.41(4) and 
t e g^ 2 /h = 0(4) Hz. However, a broadening is clearly 
seen for the |—5/2) state (where K^ 5 ^ 2 = 1.86(3) and 
t~fi 2 /h = 53(3) Hz) as well as the | —1/2) state (where 
Kq 1 ^ 2 = 1.19(5) and t~^ 2 /h = 117(6) Hz). This demon¬ 
strates spin-selective pinning, a crucial ingredient for the 
Falikov-Kimball model, the Kondo (lattice) model and 
other models for impurities or disordered systems [40, 41]. 

We now study the expansion of an interacting spin- 
mixture, by simultaneously loading atoms in the | — 9/2) 
and | — 5/2) state. During the modulation, the scattering 
length of 257(l)a 0 , where a 0 denotes the Bohr radius, 
also varies between 1280(80)ao and 217.0(l)ao [42], In 
Fig. 3c we see that the expansion is still very different for 
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FIG. 3. Spin-dependent expansion dynamics for Kg'^ 2 = 
2.41(4) (a) Whilst the effective tunneling for the |— 9/2) 
atoms is suppressed, atoms in other states are still itiner¬ 
ant and can the cloud expands, (b) Second moment w x 
of the real-space density distribution of spin-polarized non¬ 
interacting | — 9/2) (blue circles), | — 5/2) (red diamonds) and 
| —1/2) (green squares) atoms, compared to their initial val¬ 
ues. (c) Expansion of a repulsively interacting mixture of 
| —9/2) and | — 5/2} atoms. Data show mean ± s.d. of 5 (a) or 
9 (b) measurements. Solid lines are linear fits. 


the two states, however the |—9/2) component now shows 
a slight broadening. This may be caused by additional 
terms in H 0 g which have been predicted to arise owing to 
the presence or modulation of interactions, in particular 
for low frequencies [22, 23, 43-47], and is an interesting 
subject for further studies. 

A question of fundamental interest when studying ef¬ 
fective Floquet Hamiltonians concerns the evolution from 
a static to a modulated Hamiltonian as well as the stabil¬ 
ity of interacting Floquet systems [21, 48, 49]. From an 
experimental point of view, excitations created by mod¬ 
ulating the system can be observed as heating and atom 
loss, both of which are detrimental. In order to study 
how heating and losses depend on the amplitude and fre¬ 
quency of the modulation, we start with 233(10) x 10 3 
atoms in a balanced mixture of |—9/2) and |—5/2) at a 
temperature of 19(1)% of the Fermi temperature Tp. We 
turn on a three-dimensional cubic Hubbard-regime opti¬ 
cal lattice with an isotropic tunnelling of either 174(9) Hz 
or 67(3) Hz within 200 ms and then linearly ramp up the 
modulation in 100 ms. After a variable waiting time r\y 
this loading procedure is reverted (see Fig. 4a). The final 
temperature after thermalization is extracted from Fermi 
fits to the momentum distribution [34] . The measure¬ 
ments were taken at a static field of 160.44(1) G, where 
the scattering length of 194.8(l)a 0 is modulated by less 
than ±7ao for the modulation parameters of Fig. 3. 

Fig. 4b shows how the ramp-induced heating depends 
on the modulation amplitude, when setting a frequency 
of 750 Hz, as above. For small modulation amplitudes, al¬ 
most no additional heating compared to the effects of the 
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FIG. 4. Heating measurements, (a) Schematic of the ramp 
protocol, (b) Temperature after the full ramp (rw = 0), heat¬ 
ing rates and lifetimes as a function of modulation amplitude. 
Grey circles (orange diamonds) indicate a lattice with 174 Hz 
(67 Hz) tunnelling rate. The line (shaded region) show the 
results (error bars) in the dipole trap without the optical lat¬ 
tice and gradient modulation, (c) Same quantities measured 
as a function of modulation frequency, for = 1.0. Zero 

frequency indicates no modulation. Data and error bars are 
mean ± s.d. of 10 measurements (upper panels), or the results 
of linear (center) or exponential (lower) fits to the tempera¬ 
ture and atom number versus rw- 


lattice-ramp itself are observed, especially for the deeper 
lattice. For larger amplitudes, some heating becomes vis¬ 
ible at this frequency, which interestingly seems to sat¬ 
urate when K 0 9 ^ 2 > 2.4 and the effective tunneling be¬ 
comes small. A heating rate is extracted from linear fits 
to the temperature as a function of % for times up to 
300ms (after that, the temperature starts to saturate). It 
shows negligible values for most experimentally relevant 
time-scales, even for strong modulation. Exponential fits 
to the atom number measured for waiting times up to 2s 
show that modulation decreases lifetimes in the lattice. 
However, even when completely localizing one species, 
very long lifetimes of several seconds are still observed, 


which correspond to values several orders of magnitude 
larger than the interaction and tunnelling times. 

Further insight into the relevant excitations of the sys¬ 
tem can be gained by studying the dependency on modu¬ 
lation frequency. The excitation of doubly occupied sites 
is predicted to have a resonance around 1.3 kHz (2.0 kHz), 
given by the on-site interaction energy, and the width 
of the lowest band is 2.1kHz (0.8 kHz) in the shallower 
(deeper) lattice. Remarkably, directly above the sum 
of these two frequencies, where such excitations should 
be strongly suppressed [50] , the modulation ramp seems 
to cause no heating at all, see Fig. 4c. For even higher 
frequencies, the ramp-induced heating, heating rate and 
lifetime rapidly worsen for the shallower optical lattice, 
but not in the deeper one. We therefore attribute this 
feature to excitations of higher bands. Although the di¬ 
rect gap of the non-interacting band structure has a value 
of 14.7 kHz (21.8 kHz) for the shallower (deeper) lattice, 
an oscillating force can also drive higher-order processes. 

In conclusion, we have demonstrated a versatile 
method for creating widely tunable state-dependent lat¬ 
tices with minimal heating and loss, which should be 
easy to implement for many existing experimental se¬ 
tups. We have studied the static and dynamic behaviour 
of fermions in spin-dependent lattices in both real- and 
momentum-space. This method makes numerous many- 
body Hamiltonians accessible for ultracold atoms [56-59] , 
including limiting cases of the spin-anisotropic Hubbard 
model such as the Falicov-Kimball or XXZ model. The 
effects of explicitly breaking SU(2) symmetry may in fact 
already become visible at the level of short-range mag¬ 
netic correlations observed in the Fermi-Hubbard model 
so far [60-62] . In addition, interesting extensions of topo¬ 
logically non-trivial Hamiltonians can be accessed, such 
as interpolating between the Haldane and Kane-Mele 
models [63, 64]. With minor modifications our scheme 
could also be used to engineer gauge fields and spin-orbit 
coupling [51-54], where circumventing spontaneous emis¬ 
sion has been identified as a major challenge for future 
work [55] . 
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SUPPLEMENTAL MATERIAL 
Preparation of the Atomic Cloud 

We start with 1.4(4) x 10 6 non-interacting spin- 
polarized fermionic 40 K atoms which are cooled to tem¬ 
peratures of about 28(2)% of the Fermi-temperature by 
sympathetic cooling with 87 Rb in a magnetic QUIC- 
trap [35] . The atoms are transferred to an optical dipole 
trap operating at 826 nm and transferred to the required 
spin-state using radio-frequency transitions. For mea¬ 
surements with interacting spin mixtures (Figs. 3c and 
4), a mixture of | — 9/2) and | — 7/2) atoms is evaporated 
in the optical dipole trap and the | — 7/2) atoms are then 
transferred to the |— 5/2) state. For measurements with 
spin-polarized non-interacting clouds (Figs. 2 and 3b), 
the depth of the optical dipole trap is simply lowered un¬ 
til the desired atom number and a narrow momentum 
distribution is reached. 

For the measurements in Fig. 2, the atoms are ad¬ 
ditionally levitated using a static magnetic gradient 
(gravity points along the ^-direction). In Fig. 2, 
the overall harmonic confinement is u x ,y,z = 27r x 
(15.6(1), 27.8(2), 54.0(5)) Hz and in Fig. 3 it is u} x>VlZ = 
2ir x (67.8(3), 60.4(4), 233.5(3)) Hz before the expan¬ 
sion starts and <jJ x ,y,z = 27r x (0,58(1), 124(2)) Hz dur¬ 
ing the expansion. In Fig. 4 it is x> x . y . z = 2i r x 
(55(3), 52(3), 112(3)) Hz in the 174(9) Hz lattice, which 
has a potential depth of 7 Er per lattice beam (E& = 
h 2 /(2?71kA 2 ) denotes the recoil energy of 40 K with mass 
m-K in a A = 1064 nm lattice) and co Xt y yZ = 2ir x 
(66(4), 61(4), 118(4)) Hz in the 67(3) Hz lattice, which has 
a depth of 11 Er. 

Gradient and Field Calibration 

All gradients are calibrated using Bloch oscillations of 
the |—9/2) and |—5/2) atoms in a static gradient. The 
dimensionless modulation amplitude Kfi is then found 
by dividing the Bloch-oscillation frequency by the mod¬ 
ulation frequency. When working at low external offset- 
fields, the changing offset-field arising from modulation 
with a single coil (which forms part of the QUIC setup) 
leads to a non-linear dependence of the gradient on the 
driving current. A static gradient which would arise 
from the non-linearity is compensated. The non-linearity 
also introduces higher harmonics in the modulated force. 
However, it does not modify the effective tunnelling by 
more than 2%. 

In Fig. 2, the levitation of the atoms leads to a different 
external field depending on which of the two spin-states 
is supported against gravity. Therefore, driving a current 
through the coil results in a different magnetic gradient 
at the position of the atoms for the two different config¬ 
urations, which is taken into account in our calibration. 


The static and time-dependent magnetic fields are cali¬ 
brated by spin-flips from the | —7/2) to the |—5/2) states 
as well as optical absorption measurements on the D2 
line. Note that the time-dependence of the magnetic field 
could be avoided by modulating the current of more than 
a single coil, if for example a modulation of the scattering 
length is not desired. The minimal change in magnetic 
field is then given by the applied gradient multiplied by 
the size of the cloud, which is typically less than 0.1 G. 

Quasimomentum Distribution and Dipole 
Oscillations 

In the experiment, the quasimomentum distribution is 
measured by ramping down the modulation amplitude in 
10 ms and turning off the lattice in 0.5 ms, slow enough 
that higher bands are not populated in this ramp but 
fast enough that the harmonic trapping potential does 
not change the quasimomentum distribution during the 
ramp. After 15 ms of ballistic expansion, an absorption 
image of the cloud is taken. 

The dipole oscillations in the harmonic trap are initi¬ 
ated by displacing the atoms by qo in quasimomentum 
space. As a result of the tight-binding dispersion, the 
equation of motion for the quasimomentum is equivalent 
to the one of a mathematical pendulum. The square of 
the oscillation frequency is proportional to the tunnelling 
energy t and can be written as v 2 = a(qo)7r 2 v 2 Tap t i / E-r, 
where v t rap is the trap frequency without lattice. The pa¬ 
rameter a(qo) describes the effect of the anharmonicity of 
the dispersion and depends on the initial displacement. 
For < 7 q —> 0, the time evolution can be seen as an oscil¬ 
lation in a quadratic dispersion with an effective mass 
and a(qo =0) = 1. For a finite displacement, a can 
be calculated numerically and has a value of 0.89(4) for 
our case with qo = 0.31(4) 7r/a. Therefore, we expect 
a frequency of i/ 0 = 9.2(2) Hz for the lattice configura¬ 
tion in Fig. 2 without gradient modulation, which is in 
close agreement with the experimental value of 8.4(3) Hz. 
We attribute the residual deviations from the theoretical 
value to an uncertainty in the calibration of the lattice 
depths, which however cancels when considering Ks/t- 
Since the displacement in quasimomentum space is the 
same both with and without gradient modulation, the 
factor a cancels when calculating the ratio between the 
two oscillation frequencies in Fig. 2c and (v/v o) 2 — tes/t- 

As the cloud has a finite width in quasimomentum 
space, the atoms which are closer to the edge of the Bril- 
louin zone have a longer oscillation period since a(qo) 
decreases as the displacement increases. While the dy¬ 
namics of the peak of the distribution still follows the 
simple evolution described above, the motion of the cen¬ 
ter of mass depends on the exact quasimomentum distri¬ 
bution. In particular, the center of mass oscillates with 
a lower frequency than the peak. 


For the dipole oscillation data in Fig. 2c we determine 
the oscillation frequency by fitting the peak of the quasi¬ 
momentum distribution in the lattice, summed over q y 
and q z . Within this procedure we first obtain a smoothed 
quasimomentum distribution by applying a Savitzky- 
Golay filter to the raw absorption image. In a second 
step we determine the peak position of the cloud by per¬ 
forming a center of mass evaluation for data points which 
are above a threshold of 0.85 of the maximum atomic den¬ 
sity. The oscillation frequency is subsequently extracted 
from a damped sine fit function to the peak position at 
different waiting times in the shaken lattice. We obtain 
the error on the oscillation frequency by fitting the os¬ 
cillations frequency for a threshold of 0.7 and 1.0 of the 
maximum atomic density and following the same proce¬ 
dure. 


Expansion Measurements 

For the expansion data the real-space density is mea¬ 
sured without switching off the lattice or the modulation, 
with the measurement time chosen such that it corre¬ 
sponds to a zero-crossing of the magnetic field modula¬ 
tion. We measure the different spin-states separately by 
making use of the differential shift of the optical transi¬ 
tion frequency used for imaging at strong external mag¬ 
netic fields. The width of the cloud is determined by a 
Gaussian fit to the in-situ density profile. 


Heating Measurements 

We have observed that losses become very large when 
the magnetic offset-field (which, as mentioned above, also 
varies during the modulation) reaches the Feshbach res¬ 
onance at 224.21(5) G [36]. Therefore, all heating mea¬ 
surements were taken at a static field of 160.44(1) G, 
where the scattering length of 194.8(l)ao is modulated 
by less than ±7ao- After the lattice has been completely 
ramped down, we move to an offset-field of 214.8(1) G 
within 100 ms and wait for 200 ms in this configuration, 
where the scattering length is 315(2) Bohr radii, to ensure 
that any residual excitations of the cloud have thermal- 
ized. 

The heating of the atoms caused by ramping into the 
lattice and reverting the loading procedure again gives 
an estimate of the actual temperature in the lattice. In 
previous work it has been found that the heating caused 
by the first half of the ramp is lower than that of the 
reversed ramp [50, 61]. 

When measuring the heating for long waiting-times 
Tyv, we observe a saturation of T/Tp. For the data in 
Fig. 4 we fit linear slopes to times up to 300 ms, which are 
the most relevant time-scales for the majority of optical 
lattice experiments. Fitting to longer times would result 
in lower heating-rates, but the fitting function would not 
capture the time-dependence well. The 1/e- lifetimes are 
extracted from exponential fits to the atom numbers for 
waiting-times up to 2 s. 


